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Abstract 

We calculate nuclear shadowing in lepton-deuteron deep inelastic scattering, which arises from the 
double scattering of the virtual photon from both nucleons in the deuteron. The total correction to 
the deuteron structure function is found to be ~ 1% at small x, but dependent on the model deuteron 
wavefunction. The resulting increase in the corrected neutron structure function is ~ 1 — 2% for 
x ~ 0.004, which leads to a 4 — 10% decrease in the value of the Gottfried sum obtained recently by 
the New Muon Collaboration. 
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1 Introduction 



The quark structure of the nucleon is one of the most fundamental aspects of hadron physics. Deep 
inelastic scattering (DIS) of leptons from hydrogen has yielded a wealth of information on the deep 
inelastic structure of the proton. However, the absence of free neutron targets has forced one to use 
deuterium in order to extract data on the neutron structure functions. Traditionally in DIS on the 
deuteron, in which the proton and neutron are held together very weakly, nuclear effects have been 
ignored, and the total lepton-deuteron cross section assumed to be the sum of the lepton-proton and 
lepton-neutron cross sections. It is the deviation from this simple relation in the region of small 
Bjorken x (x ~ 0.1) which is known as shadowing. 

Experimentally, a deviation from linearity has been observed ||] in the so-called nuclear EMC effect 
for the ratio of DIS cross sections (or structure functions) for lepton scattering from a nucleus and 
from deuterium. A dramatic decrease in the nuclear structure function per nucleon in the region 
of small x confirmed earlier predictions Q that shadowing should be present in DIS. Furthermore, 
the shadowing was found to be only weakly dependent on Q 2 . The extraction of information about 
the difference between nuclear structure functions and those for the free nucleon from the observed 
nucleus/deuterium ratios is sensitive to any nuclear effects in the deuteron. Conclusions made about 
nucleon parton distributions based on the nuclear/deuteron structure function ratios (eg. for the 
proton antiquark distributions in the Drell-Yan process ||) at small x may have to be modified once 
shadowing is taken into account. 

A precise knowledge of the neutron structure function, F% n , is essential for the determination of the 
Gottfried sum rule, and the corresponding resolution of the question of flavour symmetry violation 
in the proton sea. It is necessary therefore to check for nuclear shadowing effects in deuterium and 
include this correction in the extraction of F% n from the deuteron structure function, Fid- Some recent 
estimates J3|, [|] have suggested a significant amount of shadowing in deuterium (up to 4%) for x ~ 0.1. 
Other calculations || have predicted a less dramatic effect 2%). 

The cross section for lepton-deuteron DIS, fig.l, is related to the forward -f*D scattering amplitude. 
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In the impulse approximation, fig. 2, the virtual photon interacts with one of the nucleons in the nucleus. 
The double scattering diagram, fig. 3, in which both nucleons participate in the interaction, is the origin 
of the shadowing in a nucleus. 



2 Vector Meson Dominance 

Hadron-Deuteron Glauber Scattering 

Glauber theory |7|, ^] for hadron-deuteron scattering gives the total hD cross section as a sum of 
the hN cross sections, and a screening term arising from the double scattering of both nucleons: 

(JhD = 2PhN + 5cThD (1) 

where 

5 °hD = -^f / d 2 k T S D (k 2 ) 

dk k S D (k 2 ), (2) 



a 



hN 



4lT 

with k = |k|. In deriving 5ahD, the assumption is made that the hadron-nucleon scattering amplitude, 
ThNi is primarily imaginary, KeJ^hN ^ Im-^/iiVj and approximately independent of k 2 for small k 2 . 
(Contributions to SohD from large k 2 will be suppressed by the deuteron form factor ^(k 2 ).) Then 
from the forward double scattering amplitude || 

SFhD = Jd 2 k T S D (k 2 ) ^ p (k 2 ) ^ n (k 2 ), (3) 

where q is the momentum of the projectile, eqn.(|2|) follows via the optical theorem: 

a = -, — r lm>. 

H 

7*D Scattering 

Assuming that the Glauber formalism can be applied to j*D scattering, the shadowing correction to 
the <y*D cross section was originally calculated in terms of the vector meson dominance (VMD) model, 
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where the virtual photon dissociates into its hadronic components (vector mesons) before interacting 
with the nucleon - see fig. 4. In this model the shadowing cross section is given by pi 



where v = p° ,u>, (j), and the photon-vector meson coupling constants are [11] 

f v a 2 M v 



47r 3 T v ^ e+e - 



(5) 



(equal to 2.28, 26.14, 14.91 for p , uj and respectively []). Writing (|j) in terms of the deuteron 
structure function, F2D Q> we have 



where now 



&°vD = / d 2 k T S D (k 2 ). (7) 

The total vector meson-nucleon cross sections, ct v n, are related to the total ttN and KN cross sections 
via the quark model, and are set to 24 mb for v = p° and u, and 14.5 mb for v = <ft (see [10, 12]). The 
deuteron form factor 5^)(k 2 ) is given by the electric monopole body form factor |l3| 



S D (k 2 ) = J™ dr(u 2 (r)+w 2 (r))j (kr), (8) 



1 Note that the fine structure constant evaluated at Q 2 — C(lGeV 2 ) is a — e 2 /A-k w 1/130, although the error 

introduced by this is probably less than that associated with using f 2 , which is obtained from the decay of meson v with 

time-like Q 2 , for the coupling to a photon with space- like Q 2 . 

2 In terms of the total cross section for the photo- absorption of virtual photons on an unpolarised deuteron, a^'n, the 



deuteron structure function is 



w K Q 2 

W2D = . o „ , 7°VD 

47T 2 Q Q 2 + V 2 



where K — \J v 2 + Q 2 is the flux of incoming virtual photons (in the Gilman convention) , so that in the Bjorken limit 



Q 2 

F2D = . , O-i'D- 
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where u(r),w(r) are the S,D-wave deuteron wavefunctions, normalised s.t. / dr (u 2 (r) + w 2 (r)) = 1, 
and where jo is the spherical Bessel function. The square of the 3-momentum transfer to the interacting 
nucleon is k 2 = + k 2 L , where k\ = m 2 N x 2 (l + M 2 /Q 2 ) 2 , and x = Q 2 /2p ■ q. 

From eqn.(||) it can be seen that the VMD shadowing correction to the deuteron structure function 
decreases as 1/Q 2 for Q 2 — > oo. 

At Q 2 = 4GeV 2 the VMD model shadowing predictions are given in fig. 5 for deuteron form factors 
obtained from several different NN potential models. By far the largest contribution (~ 80%) to 
§( V )F2D comes from the p° meson. The magnitude of 6^ F2d(%) decreases with x because the lower 
limit of the fc-integration in eqn.(Q), k m i n = ki, is an increasing function of x, and the integrand 
peaks at small values of k (~ 0.7 fm _1 ). The model dependence arises from differences in the large- 
k (~ 2 fm _1 ) behaviour of the form factor, fig. 6, which itself is determined by the small-r behaviour 
of u(r),w(r). All of the deuteron wavefunctions obtained from realistic NN potential models (namely 
Paris @, Bonn (OBEPQ) |1| and Bochum @) produce a trough in kSn(k 2 ) at k ^ 3.5 fin -1 
(because the Bessel function is negative at large kr), and a rapid fall-off with k for k ~ 6 fm _1 . Also 



shown is the model of Franco and Varma [17], which was used in [|], ||, for which the form factor, 
parameterised by a sum of Gaussians, has no laige-k tail at all. The form factor with the Paris 
wavefunction, which has the 'deepest' trough, leads to 8^ V ^F2D which is ~ 25% smaller for x ~ 0.01 
than with the Franco and Varma form factor. The trough is also responsible for the antishadowing in 
the region x ~ 0.2. 



3 Diffract ive Scattering from Partons 

At low Q 2 , it is most natural to evaluate the j*D shadowing in terms of the VMD model. At 
higher energies a parton picture may be more relevant. An alternative description of the double 
interaction mechanism in fig. 3 in the high energy limit is in terms of Pomeron (V) exchange, fig. 7. 
If the momentum transfer between the photon and nucleon is small, the nucleon will most likely 
remain intact, in which case there will only be exchange of vacuum quantum numbers. Although 
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there is as yet no QCD-based derivation of the properties of the reactions described by Pomeron 



exchange (eg. constant hadronic cross sections), there have been suggestions [18, 1£, |2(| that the 



Pomeron represents a system of gluons. (In ref.[18| hadron-hadron scattering is modelled in terms of 



gluon exchange between MIT bags, while in ref.|20| gluon-ladder techniques are used to calculate deep 
inelastic structure functions of hadrons at low x.) 

In fig. 7 the virtual photon probes the parton structure of the Pomeron, which is parameterised by 



the Pomeron structure function F 2 p [21, (defined in terms of the cross section for 7*-Pomeron 
diffractive scattering) : 

Q 2 

F ^{x) = 4^ a 7*7>- ( 9 ) 

The contribution to the F 2 d structure function from multiple diffractive scattering with V exchange 
can be written as a convolution of an exchange-P distribution function, f-p(y), with the V structure 
function: 



5 {V) F 2D {x) = f 2 dy f v (y) F 2V {x v ) (10) 

Vmin 



where 



My) = -£f \ I d2kT SD{k2) 

is expressed as a function of the momentum fraction of the nucleon carried by the Pomeron, y = 
k • q/p • q = x(l + M\IQ 2 ) M\js (M\ = p^,s = (p + q) 2 ), and we define x-p = x/y. Fig.8 
illustrates the y-dependence of fp(y), including the 1/y divergence for y — > 0. The rapid fall off with 
y is testament to the very small contribution coming from the large-y region. 

In formulating a complete description of shadowing which includes more than one mechanism care 
must be taken to avoid possible double counting. Because of this concern some authors || have 
restricted the Pomeron exchange process to the region of M\ above the highest mass of the vector 
mesons contributing to the VMD process: M\ > M^ o ~ 1.5GeV 2 , and consequently have taken the 
lower bound on the integral in eqn. (11) to be ymin = x(l + M\JQ 2 ). The VMD contribution, which 
is essentially a higher twist (1/Q 2 ) effect, may compete with that part of the diagram in fig. 7 which 
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contains low-Mx single particle intermediate states. By keeping only the leading twist piece of the 
structure function Fzp-> we can exclude this contribution since it involves extra factors of 1/Q 2 from 
the electromagnetic form factors. Nevertheless, we have tested the sensitivity of our numerical results 
to the cut-off procedure by varying M Xq from to 2 GeV 2 . For low x we find a difference over this 
range of only some 5% of the total V exchange contribution to Fid ■ For larger Q 2 the separation into 
separate Mx regions becomes irrelevant since y m in — * % m the Bjorken limit. 

For the Pomeron structure function we include contributions from the quark-antiquark box diagram, 



fig. 9a, and from the triple Pomeron interaction, fig. 9b (see refs.[23, 24|): 

F 2P {x v ) = F^\x P )+Fg?\x P ) (12) 

normalised such that 



IGthA d 2 F? ff 



dt dy 



(13) 

i=0 



where t ~ — k^, and F^ 33 is the diffractive structure function, describing semi-inclusive diffractive 
lepton-nucleon DIS, in which the recoil nucleon and the hadronic state X are separated by a large 
rapidity 



The Pomeron structure function arising from the quark box diagram, F^p X \ has been calculated 



by Donnachie and Landshoff [22 



r?{box)( \ ( 12S g 2 N sea) 0Q , , 

F 2 v p '(x P ) = 2— x v {l-x v ). (14) 



32 _ q An.*.\r-2 



The quark — Pomeron coupling constant is 1 = 3.4GeV~ [25], and we assume the same strength for 
u, d quark and antiquark — Pomeron couplings, but a weaker coupling to the strange quarks: T, q 2 = 
(10 + 2A s )/9 with X s ~ 0.5. According to the Particle Data Group pT| , the proton-proton total cross 
section a pp is approximately 40 mb. The parameter N sea is determined by the x — > behaviour of the 
nucleon sea distribution, xq sea (x — » 0) — > N sea x a . Recent parameterisations of world DIS, Drell-Yan 
and prompt photon data (2(1 12^, E8| give N sea ~ 0.15, and a approximately 0. Note that the overall 



normalisation of the r.h.s. of eqn.(14) is slightly smaller than in || due to our smaller sea parameter 



N sea (cf. N sea = 0.17 in B) and suppression of strange — Pomeron couplings. More recently, Nikolaev 
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and Zakharov [24] have calculated the box diagram contribution to 1*27? > based on a perturbative QCD 
analysis of qq fluctuations of the virtual photon. The xp dependence of their F^p parameterisation 
is the same as that in eqn.([L4|): M X /(Q 2 + M x ) 3 (since Q 2 + M x = Q 2 /x-p from the definition of x-p), 



providing the same normalisation is used (the normalisations in [22] and 2£] differ by an overall 
factor 1 — x-p). 

The triple Pomeron part of the V structure function, 

V d 2 a hp ^ hX 



7 (sp) 

2P 



(X-p) 



16vr 



a hp dtdy 



t=0J 



Fi e N a (x v ,Q 2 



(15) 



follows from 



1 d 2 F: 



2 zpdiff 



r 2N 



dtdy 



1 d 2 ahp^hx 



t=o 



a h P dtdy 



(16) 



t=o 



and the Regge theory expression for the diffractive differential cross section 

d 2 a hp ^ hX Phvjt) (3p V {t) gspjt) 1 _ 2ap ( t ) 



-y 



(17) 



dtdy 167T 

where ap(t) ~ 1 + 0.25t. In the Regge model the total hp cross section is also given in terms of the 
hadron — Pomeron couplings, f3hp- (Jhp = Phv(0)(3 p p(0). It is then evident that the combination 



1 d 2 (Jhp^hX 



cr hp dtdy 



t=0 



/V(0) 93V (0) 
I6iry 



(18) 



is independent of hadron h. From experiments on the diffractive dissociation of 7r =l= ,i^ ± ,p and p 
on hydrogen, the triple Pomeron coupling constant was found to be #3^(0) ~ 0.364 mb 1 / 2 [ 31 1 , 
independent of t, and indeed of the hadron type h. 



For the sea part of the nucleon structure function, F^ 



hx(u s + u + d s + d+2(s + s)/5)/18, we use 



recent parameterisations of the data at Q 2 = 4GeV 2 [27, 25]. In the calculation of ref.||, a constant 
value of 0.3 was used for F^ff together with an empirical low-Q 2 dependence p2]] . With the above 
triple Pomeron coupling constant, eqn.([l5|) gives a 3V component which is about 40% smaller than 
that obtained in However, this is not very significant for the total Pomeron structure function, 
since F^p is very much smaller than the quark-antiquark 'box' contribution, F^p° x \ fig. 10. 
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The scaling behaviour of the "P-exchange mechanism is determined by the scaling behaviour of the 
V structure function, and from eqns.(p^) — (18) it is clear that will scale as Q 2 — > oo. At Q 2 = 



4GeV 2 , fig. 11a shows the individual 'box' and 3V contributions to S^F2D, with the deuteron form 
factor obtained from the Bochum wavefunction. The dependence of 8^ F2D on S£>(k 2 ) is illustrated in 
fig. lib. Again, as in the case of the VMD model, the large-/c negative tail of the form factor produces 
a large (some 30-40%) difference between different models for x ~ 0.05. For x ~ 0.2 the presence or 
absence of antishadowing will be determined by the model deuteron wavefunction. 

4 Shadowing by Mesons 

Another potential source of shadowing arising from the double scattering mechanism is one which 



involves the exchange of mesons, fig. 12. It has previously been suggested []32j that this leads to 
substantial antishadowing corrections to F2d(x). The total contribution to the deuteron structure 
function from meson exchange is written 

5 {M) F 2D (x) = £ / d v fM (19) 
ft Jx 

where \x = n,p,u,a and y = k ■ q/p ■ q = (ko + ki)/rriN and x^ = x/y. For the virtual meson 
structure function, i*2u, we take the parameterisation of the (real) pion structure function from Drell- 
Yan production p3[ |. The exchange- meson distribution functions /^(y) are obtained from the non- 
relativistic reduction of the nucleon — meson interaction: 

The deuteron wavefunction is defined by 

tf(p,J.) = (u(p) - w(p)^^j x(% (21) 



where u(p) and w(p) are its S and D-w&ve components, normalised so that J dp p 2 (u 2 (p)+w 2 (p)) = 1, 
with p = p/p and p = |p|, and 5i2 is the tensor operator: 5i2(p) = 3 a\ ■ p (X\ ■ p — a\ ■ 02- The 
deuteron spin wavefunction is denoted by x[ z > where J z is the total angular momentum projection. 
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In eqn.(20), k 2 = Uq — k 2 , where fco = mp — Jrn 2 N + p 2 — Jmfj + p' 2 is the energy of the off-shell 



meson, and k = p — p' is its 3-momentum. 

The nucleon — meson interactions are given by JlS] 



VnNN 
VpNN 



I, 



ttNN 

ml 



o\ ■ k a 2 • k 



(22) 



2 

9pNN 



1 + 

9pNN fpNN 



3q 2 



2777^ 8777^ 171 ° 2 4m 2 N 



k 2 o"i • k o"2 • k 



4m 2 N 



2rriN 



k 2 k 2 (Ji • k (72 • k 
<7i • a 2 1 



m N 



m N 



m N 



(23) 



+ 



VujNN 
VrrNN 



f 2 

JpNN 

4m 2 N 



2 

9u)NN 



-0\ '(72 k + 0"! • k CT2 • k 



1 + 



2 

~9 a NN 



3q 2 k 2 

2777^ 8777^ 

q 2 k 2 

+ 



4777^ 



+ 



0"! • k 0"2 • k 

4777^ 



2777^ 8777^ 



(24) 
(25) 



where q = ^(p + p')- Terms proportional to S • k x q, where S = o\ + a 2 , are omitted as they do not 
contribute to fp(y)- 

Evaluation of eqn,(|2G|) requires the identities: 

\^2&(p,j z ) *(p',j z ) = ^ [u[p) u{ P ') + w{p) w(p') p 2 (co S e) p 2 (co S e')] 



+ (j) dependent terms 

= \Y,^ip,Jz)(Tl-(T 2 ^(p',J z ) 



(26) 



Jz 



4vr I 3 

1 
71 



k 2 - 2 p p' sin# sin#' u(p) 
4 p p' cos 6> cos 9' sin 2 6>' + 4 p' 2 cos 2 9' sin 2 6>' 



o (P 2 + P ) p 2(cos 0') + 2 (p 2 cos 2 9 + p' cos 2 0') P 2 (cos 

o 



+ - p p' cos 9 cos 9' P 2 (cos9') 
3 



71 (p) 



1 
71 



/ / 2 2 2 2 

4 p p cos 9 cos 9 sin 9 + 4 p cos sin 
(p 2 + p 2 ) P 2 (cos 9) + 2 (p' 2 cos 2 0' + p 2 cos 2 0) P 2 (cos 

ro(p) u(p') 



+ - p' p cos 9' cos 9 P 2 (cos 
3 
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- (p co&e+p' cosO') 2 P 2 {cos6) P 2 COS V 
3 L 

- 2 (p 2 sin 2 9 + p' 2 sin 2 9') P 2 (cos9) P 2 {cos9') 
+ 3 (p 2 cos 2 6> sin 2 6 + pp' cos 9 cos 9' sin 2 9) P 2 (cos9') 
+ 3 (p /2 cos 2 sin 2 9' + p' p cos 9' cos 9 sin 2 9') P 2 (cos9) 

( w(p) w(p')\ (27) 
+ dependent terms. 
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+ — p p' cos 9 cos 9' sin 2 9 sin 2 9' 



The terms in eqns.(|26|)-(j27|) which depend on the azimuthal angle (4>) vanish after integration. The 
factors cm are due to isospin: c w = c p = 3,c w = c a = —1. The fiNN vertex form factors T^nn^ 2 ) 
are parameterised by a dipole form 

/A 2 - m 2 \ 2 

^NN(k 2 ) = [-^Z^t) > ( 28 ) 



with the high-momentum cut-offs ranging from ~ lGeV in models with soft form factors [34, 16| 
to ~ 1.7 — 2GeV when hard form factors are employed fl~5|| . 

Fig. 13 shows the individual meson exchange contributions to 5 <yM ^ > F 2 £ ) , for the wavefunction of 
the Bonn model, and with a universal dipole cut-off of A^ = 1.7GeV. As could be expected, pion 
exchange is the dominant process. We also include the fictitious a meson, but with a mass (~ 800 
MeV) that is larger than that used to represent 2ir exchange in NN scattering. Both of these produce 
antishadowing for small x. The exchange of vector mesons (p, oS) cancels some of this antishadowing, 
although the magnitude of these contributions is smaller. In fact, for A^ ~ 1.3 GeV all contributions 
other than that of the pion are totally negligible. 

Fig. 14 shows the dependence of the total 5^ M ^F 2 o on A M for the Bonn model wavefunction. There 
is approximately a factor of 2 difference between the amount of shadowing with soft (A^ = lGeV, 
lower dotted line) and hard (A^ = 1.7GeV, upper dotted line) form factors. In lepton-nucleon DIS it 
is well known p5| that the meson cloud of the nucleon, with a hard fj,NN form factor, gives nucleon 



sea distributions that are several times larger than the empirical ones. In fact, to be consistent 
with the lepton-nucleon DIS data A M must be ~ 0.8 — 0.9GeV. We also consider the effect of the 
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model momentum-space deuteron wavefunction on 5^ M ^ F 2 £> ■ Although the model wavefunctions differ 
substantially at large momenta (p ~ 2 fm -1 ), this variation will be largely suppressed by the /j,NN 
form factor. The Bochum and Paris wavefunctions are generally larger than the Bonn wavefunction, 
and this is reflected in a larger 5^F 2D . 

We also comment here on the issue raised in the previous section, namely double counting, this 
time between the meson exchange and the other mechanisms. It should be clear that since the V 
contribution involves the exchange of vacuum quantum numbers, there will be no interference between 
this and the exchange of pseudoscalar pions or vector mesons. The scalar a meson, introduced as an 
effective description of two-pion NA excitations, does not correspond to actual exchange of a spin 
particle. By restricting the meson structure function to only the leading twist component (our F 2 m 
is determined at Q 2 = 25 GeV 2 where this assumption is reasonable) we may view the VMD process 
as a description of higher twist effects. Still, imposing any low-Mx cut on the meson exchange term 
has numerically insignificant consequences, largely because F 2fJL (x/y) — ► as y — > x. 

5 Combined Shadowing Effects and the Gottfried Sum Rule 

The total deuteron structure function is defined by 

F 2D (x) = F 2p (x) + F 2n (x) + 5F 2D (x), (29) 

where the shadowing correction is a sum of the VMD, Pomeron and meson exchange contributions: 

5F 2D (x) = 8^F 2D {x)+5^F 2D (x)+5^F 2D {x). (30) 

In fig.15 we compare the contributions to 5F 2 £>(x) from the three mechanisms considered. For x ~ 
0.1 the magnitude of the (negative) Pomeron/VMD shadowing is larger than the (positive) meson- 
exchange contribution, so that the total 5F 2 n is negative. The fact that shadowing is present in this 
region of x does not depend on the model deuteron wavefunction. For larger x (« 0.1 — 0.2) there is 
a small amount of antishadowing, which is due mainly to the VMD contribution. The dependence of 
the total shadowing correction on the deuteron wavefunction and on the (iNN form factor is shown 
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in fig. 16 for Q 2 = 4 GeV 2 . We point out that the magnitude of 5F 2 d{x) is about 4 times smaller than 
that obtained in reference and about 2 times smaller compared with the result of reference ||. 
The most important reasons for our smaller results are the inclusion of meson exchange contributions 
which produce antishadowing at small x, and the use of realistic deuteron wavefunctions which lead 
to smaller V exchange and VMD contributions. 

Recently the New Muon Collaboration (NMC) has measured F 2p and F2D [36, 37 1 down to very 



small values of x (= x m i n = 0.004). The neutron structure function was then extracted from Fid in- 
order to test the Gottfried sum rule pgj ]. However, by assuming that 

F 2D (1 - (F 2D /F 2p - 1)) _ 

(1 + (F 2D /F 2p - 1)) = 2F *~ F ™ = ^v-F2n)NMC (31) 

the NMC ignored any nuclear shadowing effects in D which may alter the F 2n values. The actual 
difference between the p and n structure functions should be 

F 2p -F 2n = (F 2p -F 2n ) NMC + SF 2D , (32) 

and this is shown in fig. 17. The dotted line is a best fit to the NMC data, and includes the small- x 
extrapolation used in [37|: 



F 2p (x) - F 2n (x) ^5 ax* (33) 

with a = 0.21,(3 = 0.62. The other curves include the shadowing corrections to the NMC data 
parameterisation. It is not clear whether F 2p — F 2n will become negative at x ~ 0.004, and it will be 
interesting to see whether this cross-over occurs when additional data at smaller x become available. 
The Gottfried integral 

So(M) = ftf S^LW (34) 



d(logx') (F 2p (x') - F 2n (x')) 



is given in fig. 18 for x down to 0.004. In the naive quark model, Sg(0,1) = 1/3. Ignoring nuclear 
effects, the NMC obtained SG(x m inA) = 0.229. From the unmeasured region (x < 0.004), using the 
above extrapolation, the contribution was found to be Sq(§, x m i n ) = (a/P) x P min = 0.011. With the 
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Model 


a (3 




$G\-Emini 1) 


S G (0,l) 


NMC @ 


0.21 0.62 
0.109 0.5 


0.011 
0.014 


0.229 


0.240 ± 0.016 
0.243 


Bochum (A„ = 1.3GeV) 


0.043 0.5 


0.005 


0.222 


0.227 


Paris (A M = 1.3GeV) 


0.052 0.5 


0.007 


0.224 


0.230 


Bonn (A^ = 1.3GeV) 


0.011 0.5 


0.001 


0.215 


0.217 


Bonn (A„ = l.OGeV) 


0.002 0.5 


0.000 


0.214 


0.214 


Bonn (A^ = 1.7GeV) 


0.019 0.5 


0.002 


0.217 


0.219 



Table 1: Small- x extrapolation parameters for i*2p — -?2n( = cxx") and the contributions to the Gottfried 
sum from different x-regions. 

conventional Regge theory assumption that /3 = 0.5, Sa(0,x m i n ) would be 0.014. In table 1 we give 
the values of Sg including shadowing corrections, and also the x < x m i n extrapolation parameters. 
For simplicity we take (3 = 0.5, and adjust a to achieve a smooth transition between the x > x m i n 
and x < x m in regions. The overall correction to the NMC value for Sq(0, 1) is found to be between 
-0.010 and -0.026. This is to be compared with -0.07 to -0.088 obtained in g |, g|]. 

As a fraction of the total F2d{x) |37jl , the shadowing correction amounts to (0.5-1.0%, 0.4-0.8%, 
0.0-0.3%) at x = (0.004,0.01,0.1), while the antishadowing is less than 0.2% of F 2D at x » 0.2. 

In fig. 19 we show the ratio of neutron structure functions with and without shadowing corrections, 

F 2n = 1 _ SF2D ( l + (F 2 n/F 2p ) NM c \ (35) 
{F 2n )NMC F 2 n y {F 2n / F 2p )nmc ) 

where the NMC neutron/proton ratio was defined as {F 2n j F 2p )nmc = F 2 d/F 2p — 1. There is an 
overall 1 — 2% increase in the neutron structure function due to shadowing for x ~ 0.01. 

Finally, we illustrate in fig.20 the dependence upon Q 2 of the total shadowing correction, 5F 2 d(x, Q 2 ). 
As expected, the VMD term vanishes rapidly with increasing Q 2 , leaving the two scaling contributions 
from V and meson exchange to largely cancel each other for Q 2 ~ 25 GeV 2 . However, we should add 
a note of caution about comparing shadowing corrections at very large values of Q 2 . In the parton 
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recombination model g, [3^, 40] the fusion of quarks and gluons from different nucleons introduces 
additional terms in the Altarelli-Parisi equations governing the QCD evolution of the parton 
distributions. At very small x and large Q 2 , such as those attainable at HERA energies, this can 
lead to significant corrections || to the 5F2d(x,Q 2 ) evolved without these terms, although the exact 
magnitude of these is sensitive to the small- x behaviour of the input nucleon gluon distribution. For 
the moderate range of Q 2 and not too low x values in fig. 20, however, we expect the indicated Q 2 
behaviour to be reliable. 



6 Conclusion 

In summary, we have estimated the nuclear shadowing in lepton-deuteron DIS from the double scatter- 
ing mechanism in fig. 3. Our approach is similar to that of refs.[|| and ||, in describing the interaction 
in terms of the VMD model, together with Pomeron (V) exchange at larger Mx- However we have 
also included contributions from the exchange of mesons which effectively cancel as much as half of the 
shadowing from the VMD/'P-exchange mechanisms alone. Numerically, there is some dependence on 
the model deuteron wavefunction, and also on the meson-nucleon form factor for the meson-exchange 
process. The net effect is a ~ 1% reduction of F2D for x ~ 0.004, or equivalently a ~ 2% increase 
in the neutron structure function over the uncorrected Fi n - Consequently, the shadowing correction 
to the Gottfried sum S G (0, 1) is between -0.010 and -0.026 (or about 4 and 10% of the NMC value), 
which is about 5 times smaller than in previous estimates. 

To accurately test the descriptions of shadowing in the deuteron it is necessary to obtain model- 
independent information on the neutron structure function at low x. Even at HERA energies this is 
not possible with electron scattering alone. However, when combined with high-precision data from 
neutrino-proton experiments the individual flavour distributions can be determined, and the neutron 
structure function inferred from charge symmetry. For this to happen, however, the statistics on the 
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neutrino data need to be improved, and the range extended to smaller x. 
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Figure captions 



1. Lepton-deuteron deep inelastic scattering. 

2. Virtual photon-deuteron scattering in the impulse approximation. 

3. Double scattering of virtual photons from the deuteron, in which both nucleons take part in the 
interaction. 

4. Double scattering mechanism in the vector meson dominance model. The virtual photon 
dissociates into a vector meson which then scatters from the nucleon. 

5. VMD contribution to the total deuteron structure function at Q 2 = 4 GeV 2 , for different model 
deuteron form factors. 

6. Deuteron form factor, as defined by eqn.(||). 

7. Pomeron exchange contribution to the double scattering mechanism (the Pomeron is denoted 
by the zig-zag). 

8. Exchange- Pomeron distribution function, \f-p(y)\. 

9(a). Quark-antiquark box contribution to the Pomeron structure function, where the Pomeron 
couples to the virtual photon via a quark-antiquark pair. 

(b). Triple- Pomeron contribution to the Pomeron structure function. 

10. x-p dependence of the 'box' and 3"P contributions to F%p(x-p) for the quark distribution function 
parameterisations of Owens |27j and Morfin and Tung |28| at Q 2 = 4 GeV 2 . 

11(a). Quark-antiquark box and 3V contributions to the total deuteron structure function. The 
deuteron form factor is given by the Bochum model wavefunction. 

(b). Deuteron form factor dependence of the Pomeron exchange contribution to the deuteron 
structure function. 

12. Double scattering mechanism with meson exchange. The dotted line represents mesons 
7T, p, iv, a. 

13. Individual meson exchange contributions to the deuteron structure function, for the wavefunc- 
tion of the Bonn (OBEPQ) model with a universal form factor cut-off = 1.7 GeV. Note the mass 
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of the effective a meson is « 800 MeV. 

14. Deuteron wavefunction and fiNN form factor dependence of the total meson exchange cor- 
rection. The Bochum (solid) and Paris (dashed) curves are evaluated with A M = 1.3 GeV, while the 
Bonn (dotted) curves have A M = 1.0, 1.3 and 1.7 GeV, with the larger cut-off giving more overall 
antishadowing. 

15. Comparison between the VMD, Pomeron and meson exchange corrections to the deuteron 
structure function at Q 2 = 4 GeV 2 (for the Bochum wavefunction, and a form factor cut-off = 1.3 
GeV for the meson exchange process). 

16. Deuteron wavefunction and /iNN form factor dependence of the total shadowing correction. 
For the Bochum and Paris curves A^ = 1.3 GeV, while the Bonn curves are calculated with A^ = 
1.0, 1.3 and 1.7 GeV, with the larger cut-off giving less overall shadowing. 

17. Difference between the proton and neutron structure functions, with shadowing corrections to 
the NMC data at Q 2 = 4 GeV 2 . 

18. Gottfried sum, with shadowing corrections to the NMC data. 

19. Neutron structure function ratio, with and without shadowing corrections (at Q 2 = 4 GeV 2 ). 

20. Q 2 dependence of the total shadowing correction to F2D- Curves represent the shadowing 
corrections at 4, 10 and 25 GeV 2 , with the Bochum model wavefunction and with A = 1.3 GeV. Also 
shown is the correction at Q 2 = 25 GeV 2 without the VMD contribution. 
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